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Abstract. The evolution of a solar model is initialized 
with homogeneous models of either, pre-main sequence (P- 
models) or, zero-age main sequence (Z- models). The zero- 
age of a solar model is conventionally referenced as the 
time where the nuclear reactions just begin to dominate 
gravitation as the primary source of energy. Fixing the 
physics, we found that the structure of P- and Z-models 
computed with the same physics are almost similar soon 
after the exhaustion of their convective core. This simi- 
larity gives a connection between the age of the Sun £© 
and the time t ca \ elapsed in the calculation of calibrated 
solar models. We found that a Z-model calibrated with 
t ca \ = <© and a P-model calibrated with t ca \ = £© + 25 My, 
are indistinguishable at the relative accuracy level of a few 
10" 4 . 
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1. Introduction 

The evolution sequences of a solar model are initialized, 
either with a chemically homogeneous pre-main sequence 
model (P-models) powered only by gravothermal energy, 
or with a fictitious homogeneous zero-age main sequence 
model (Z-models) powered only by thermonuclear energy. 
Fixing the physics, whatever they start from, calibrated 
solar P- or Z-models have, by definition, the same radius, 
luminosity and surface mass fraction of heavy elements. If 
one is not concerned with specific pre-main sequence phys- 
ical processes, e.g. primeval mass and angular momentum 
losses or lithium depletion, the amount of cal culati ons is 
half as large for Z-models (see e.g. Morel et al. 1997 ). The 
initial internal structure of initial P- and Z-models differs, 
despite the same surface parameters at solar age. It is not 
obvious that the internal structure of the resulting cali- 
brated models is alike and, if so, is it just for the resulting 
models or also for parts of their previous evolution? It is 
important also to establish a relationship between the age 

Send offprint requests to: P. Morel 
Correspondence to: morel@obs-nice.fr 



of the Sun and the sum of the time steps elapsed in the 
calculation. As an example, to take the time of the pre- 
main sequence evolution into account, Weiss & Schlattl 



(1998) have evolved their calibrated solar P-models 40 My 
beyond the age they have assigned to the Sun. 



A preliminary discussion is given in Morel et al. (1998). 
The paper is organized as follows: in Sec. we recall how 
the age of a solar model is defined with respect to the age 
of the Sun. Section || is devoted to our analysis. Results 
are given in Sec. || and we conclude in Sec. ||. 

2. The solar age 

In a sequence of evolving P-models, the zero-age main 
sequence (ZAMS) model is usually defined as the first 
model where e nU c, the nuclear energy generation just be- 
gin to dominate e g , the gravothermal energy liberation 
(Guenther & Demarque 1997 ). That definition influences 
the various age definitions. Here, with this definition, for 
a sequence of solar P-models the age of the calibrated P- 
model will be defined as the evolution time i.e., the sum of 
the time steps, from the ZAMS model to the present day 
model. As we shall see that definition influences slightly 
the various age definitions. 

Conventionally, the solar age £© is the time it has taken 
the Sun to evolv e from ZAMS to present day (Guenther & 
Demarque 1997 ). According to a reasonable hypothesis, i© 
needs to be consistent with the age of the oldest meteorites 
4566 ±5 My (Bahcall et al. 1995). Guenther & Demarque 
( 1997 ) have argued that the radioactive clocks of these 
meteorites are zeroing during the last high-temperature 
event in the primordial solar system nebula. This has oc- 
curred, according to Guenther (1989), 40 ± 10 My before 
the ZAMS, therefore Guenther & Demarque (loc. cit.) esti- 
mate the "meteoritic solar age" to be £© m = 4530 ±40 My. 
The determination of i© has been recently revisited using 
new methods of investigation including statistical argu- 
ments based on helioseismological data. In their pioneering 



work Weiss & Schlattl (1998) investigated how the known 



deficits of the solar models would translate into an age, 
if the best-fitting model's age is assumed to be the solar 
age. They found a "helioseismic solar age" in the interval 
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4650 My JS i Q h ^ 5650 My. Recently another analysis of 
Dziembowski et al. (1999), based only on central helioseis- 
mic data, leads to £ Q h = 4660±110My, a value marginally 
consistent with i0 m . 



3. Solar P- and Z-models 

Hereafter we shall note i e v the evolution time of any 
model. For a calibrated solar model we shall note i ca i = 
t cv . We shall extend these notations to the Sun itself by 
setting t cv — at the time of the last high-temperature 
event in the primordial solar system nebula, thus i ca j = 
t Qnl + 40 My (Guenter loc. cit.). 

For a sequence of P-models we conventionally set 
t cv = at the ignition of deuterium; that occurs as soon 
as the temperature reaches T « 0.5 MK. In such phys- 
ical conditions the model is fully convective and then 
homogeneous. At t cv ss 0.3 My, all the initial 2 H is al- 
ready converted into 3 He via the reaction 2 H(p, 7) 3 Hc 
and < e nuc << e g . Therefore, long before the ZAMS, the 
abundance of 3 He is enhanced by a factor of rj 4 (in mass) 
with respect to its primeval value, namely from rj 2.3 10 -5 
to « 8.6 10~ 5 . At i cv w 1.2 My, the central temperature 
approaches T ca \ ~ 3MK, and the PP burning of 1 H be- 
gins in the core. Soon after, due to the sudden energy lib- 
eration, the temperature increases, the opacity decreases 
and the core becomes radiative. At t cv rj 13 My the in- 
nermost limit of the solar convection zone has receded to 
about its present day location i?zc ~ O.7i?0; the central 
temperature has jumped to T ca \ « 13 MK. The zero-age 
main-sequence occurs at t cv = £zams ~ 25 My just when 
£nuc <5 eg- At i cv w 37 My an extra energy generation due 
to the CNO burning of 12 C into 14 N creates a convective 
core. It lasts until t cv as 90 My when the CNO reactions 
reach their equilibrium state. Beyond t cv <^ 50 My more 
than 99% of the energy generated have a nuclear origin. 
After the exhaustion of the convective core, the model 
evolves quietly without fundamental modifications of its 
structure until present day. 

Therefore, for a calibrated solar P-model r p = t ca \ 
writes: 



^ZAMS + t(i 



(1) 



For a sequence of Z-models we conventionally assign 
t cv = OMy to the first model. In such a model the energy 
only comes from nuclear reactions i.e., e g = 0. It is an ide- 
alized model, because such a homogeneous state cannot 
be issued from any previous PMS evolution. As the time 
goes on, the nuclear burning and the gravity relax the ini- 
tial inconsistencies of chemical composition, temperature, 
density and pressure. At t cv = a solar Z-model presents 
a convective core which lasts until t cv w 65 My, when the 
CNO bi-cycle begins to work at equilibrium. Then, as for 
the P-models, after the exhaustion of the convective core, 
the model changes quietly until present day. 



The definition of the age given in Sec. g does not ap- 
ply to Z-models, since the gravothermal energy generation 
never overcomes the nuclear one. Thus, the time of the 
ZAMS is undefined for a sequence of Z-models. To con- 
nect £0 and t z = i cv , we infer that a P- and a Z-model 
computed with the same physics are almost identical as 
soon as the nuclear reactions start to work at equilibrium. 
That occurs « 10 My after the exhaustion of the convec- 
tive cores. Beyond this epoch the models will remain alike. 
As illustrated in Fig. [y, calibrated P- and Z-models will 
have the same age if the same amount of time is elapsed 
between present day and the epochs r cp and r cz of convec- 
tive core exhaustion: r p — r cp = t z — r cz . Owing to Eq. (|lj), 
t, writes: 



t z = r n - r, 



cp 



*ZAMS + £© — T cp + T c , 



(2) 



4. Results 



Basically, the physic s used in the calculations is the same 
as in Morel et al. (1997), it uses OPAL opacities and 
equation of state. The microscopic diffusion of chemicals 
is allowed for. The models are calibrated, within a rel- 
ative accuracy better than 10 -4 , by adjusting: the ratio 
l/H p of the mixing-length to the pressure scale height, 
the initial mass fraction Y\ nl of helium and the initial 
mass fraction (Z/X)^ of heavy elements to hydrogen, 
so that, at present day, the models have the luminosity 
L Q = 3.846 10 33 ergs" 1 , the radius R Q = 6.9599 10 10 cm, 
and the mass fraction of heavy element to hydrogen 
(Z/X)q = 0.0245. We have taken into account the most 
important nuclear reactions of PP+CNO cycles with the 
species 2 H, 7 Li, 7 Be at equilibrium. As 2 H is set at equi- 
librium, the protosolar 2 H is included into the initial 3 He; 
one has for the initial isotopic ra tio 3 H e/ 4 He = 4.19 10~ 4 
(in number) (Gautier & Morel 1997 ). The recently up- 
dated rates of Adelberger et al. ( 1998 ) are used. The mod- 
els h ave been computed using the CESAM code (Morel 
1997 ). The evolution of a Z-model necessitates about 40 



time steps whereas 90 for a P-model. Around epochs of 
convective core exhaustion and ZAMS, the time step is 
refined in order to define these instants with an accuracy 
better than ±2 My. 

We have adjusted r p for the calibrated P-model Sp 
so that its age was t Q = 4530 My. The ZAMS occurred 
at tzAMS = 25 ± 2 My and the convective core was ex- 
hausted at r cp = 90 ± 2 My. With the calibration param- 
eters obtained for Sp, we have computed the evolution of 
the Z-model Sz. We have found that the exhaustion of its 
convective core occurred at t cz = 65 ± 2 My. According 
to our analysis, after t z = t Q ± 6 My, given by Eq. (||), 
the Z-model Sz is expected to be calibrated and to have 
a structure close to the P-model Sp. That is indeed the 
case, as seen in Table |Tj and Fig. ^, as far as the global 
parameters and sound speed profiles are concerned. 

Figure |3| shows diagrams of effective temperature and 
luminosity versus time. As expected, an abscissa shift of 
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Fig. 1. Schematic illustrations of chronologies with £ Q = t Qm for the Sun, P- and Z-models; the times are in My. 



Table 1. Comparison of global characteristics of solar P- 
models (labels p and P) and Z-models (labels z and Z). 
The subscripts " s " and " c " respectively designate a surface 
and a center value. The data are for the time t cv = 
(first group), for the time 10 My after the epoch t cv of 
convective core exhaustion (second group) and for t ca \ i.e. 
present day (last group). The units are My for time, 10 7 K 
for the temperatures and gem -3 for the densities. 



Sp 



SP 



sz 



Y a 


0.2731 


0.2731 


0.2720 


0.2720 


(Z/X) s 


0.0278 


0.0278 


0.02775 


0.02775 


tcv 


90 ±2 


65 ±2 


90 ±2 


64 ±2 


Y a 


0.2724 


0.2725 






z s 


0.0196 


0.0196 






Rzc 


0.7220 


0.7221 






T c 


1.347 


1.349 






Pc 


81.35 


81.31 
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0.2781 
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Fig. 2. Relative differences in sound speed between the 
P and Z-models Sp and Sz; dashed: 10 My after the con- 
vective core exhaustion, full: present day. 



+25 My for the Z-model Sz allows to superimpose the loci 
of Sz and Sp beyond t GV « 100 My. We also found that 



the value of t c 



25 ± 4 My is almost insensitive 



to the required solar age, as seen in Table [l] from the 
comparison of the P-model Sp (resp. Z-model Sz) of age 
4530 My with the P-modcl SP (resp. Z-model SZ) of age 
4660 My. At solar age t Q = 4530 My (resp. i = 4660 My) 
the differences between the models Sp and Sz (resp. SP 
and SZ) are below 10~ 4 i.e. the level of accuracy of the 
calibration. All these results validate our analysis. 

As a matter of comparisons, we have calibrated the so- 
lar Z-models Sg and Sw for ages £© = 4500 My (Guenther 
fc De marque [1997D and t = 5100 My (Weiss & Schlattl 
1998), respectively. Figure shows the relative difference 
between the sound speed of models and the seismic sound 
speed experimental results of Turck-Chieze et al. (1997). 



The best fit appears to be for the models SP (and SZ), of 
age i© = £oh = 4660 My, i .e. the helioseismic solar age of 
Dziembowski et al. ( |l999| ). This means that the physics 
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Fig. 3. (a) Effective temperature and (b) luminosity with 
respect to time, for the solar P- and Z-models, Sp (full) 
and Sz (dashed); (c) and (d) the same with a +25 My 
abscissa shift for Sz. 




R / R, 

Fig. 4. Relative difference in sound speed between the 
Sun and models, Sp (thin, full), Sz (thin, dashed), SP 
(heavy, full), SZ (heavy, dashed), Sg (dash-dot) and Sw 
(dotted). 

used in solar models leads to deficits which arc equivalent 
to an evolution slower than indicated by tQ m . 

5. Discussion and conclusions 

Fixing the physics we have compared calibrated solar 
models initialized with a homogeneous zero-age main se- 
quence model (Z-model) and with a pre-main sequence 
model (P-model). We have verified that the models which 
evolved from these different initial models merge into the 
same structure about 10 My after the end of the convec- 
tive core phase, if the times, i ca i, of evolution used in each 
calculation, are shifted by an amount of 25 My. That sim- 
ilarity allows us to connect t ca \ to the solar age, i.e. the 
time it has taken the Sun to evolve from the time (ZAMS), 



where nuclear reactions just begin to dominate gravita- 
tion as the primary energy source, to present day. For 
the calibrated Z-model we found t ca \ — i©, while for a 
P-model t ca \ — i© + 25 My. We emphasized the fact that 
these two relations between t ca i and t Q are only valuable 
when the zero-age reference corresponds to the time when 
the nuclear reactions just begin to be the primary energy 
source. If the epoch of the ZAMS is defined at the instant 
where 99%, instead of 50%, of energy comes from nuclear 
reactions, the shift becomes 50 My, then, for a P-model 
^cai = t Q + 50 My and i ca i = t Q + 25 My for a Z-model. 
Though 25 My represent only 0.5% of the solar age, it is 
of the same order of accuracy as present day solar models. 
The basic idea of our analysis is based on the similarity 
between the models soon after the convective core exhaus- 
tion; it is useless for stars of masses greater than about 
1.2 Mq, which exhibit a convective core on the main se- 
quence. Another analysis is needed for these stars, that 
may be of importance in asteroseismology to differentiate 
between age and t ca \ for modeling the COROT targets 
(Baglin|l998l). 
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